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Abstract 

We compute the contribution to the three-point temperature cor- 
relation function of the Cosmic Microwave Background coming from 
the non-linear evolution of Gaussian initial perturbations, as described 
by the Rees-Sciama (or integrated Sachs-Wolfe) effect. By expressing 
the collapsed three-point function in terms of multipole amplitudes, 
we are able to calculate its expectation value for any power spectrum 
and for any experimental setting on large angular scales. We also give 
an analytical expression for the rms collapsed three-point function 
arising from the cosmic variance of a Gaussian fluctuation field. In 
the case of COBE DMR, we find that the predicted signal is about 
three orders of magnitude below that expected from the cosmic vari- 
ance. 
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After the detection by COBE DMR (e.g. Smoot et al. 1992, Bennett 
et al. 1994) of anisotropics in the Cosmic Microwave Background (CMB), 
several other measurements of anisotropics at different angular scales have 
been announced. This new set of observations has provided an extra tool 
to improve our understanding of the structure formation process. The usual 
analysis in terms of the two-point correlation function of the COBE data 
fixes the rms amplitude of the primordial perturbations. A comparison with 
measurements at smaller angular scales encodes information about the mat- 
ter content of the universe and other cosmological parameters to which these 
are sensitive. Furthermore, an analysis of the three-point correlation func- 
tion of the data, as it has recently been performed for the COBE DMR first 
and two-year data (Hinshaw et al. 1994, 1995), can provide useful clues 
about the statistical properties of cosmological perturbations. Further anal- 
yses probing the statistical nature of temperature fluctuations have been 
performed on the COBE maps by Smoot et al. (1994), Kogut et al. (1995) 
and Torres et al. (1995). These studies are of particular relevance as they 
could help to distinguish among the two presently preferred theories for the 
origin of primordial fluctuations: inflation and topological defects. 

Inflationary models predict a quasi-Gaussian distribution of density per- 
turbations. This is due to the fact that they are originated by the quantum 
fluctuations of a very weakly coupled scalar field, the inflaton. The effect of 
the small non-linearities in the inflaton dynamics and the mean value of the 
resulting imprints on the three-point function of the gravitational potential 
(see Falk, Rangarajan & Srednicki 1993, Gangui et al. 1994, Gangui 1994) 
are several orders of magnitude smaller than the typical values expected 
for a particular realization of an ensemble of Gaussian universes, the so- 
called cosmic variance (Scaramella & Vittorio 1991, Srednicki 1993). As the 
anisotropics observed at large angular scales are essentially determined by the 
fluctuations in the gravitational potential, the CMB three-point correlation 
function at the scales proved by COBE is also expected to be Gaussian. An 
analysis of the rms skewness of inflationary models giving rise to isocurvature 
baryon perturbations, by Yamamoto and Sasaki (1994), yields values smaller 
than those obtained for curvature perturbation models. On the other hand, 
topological defects are the typical example of non-Gaussian distributed per- 
turbations. However, for many observations, the relevant object is not the 
individual effect of each particular defect, but the superposition of many 
of them, which results in a nearly Gaussian pattern (e.g. Gott et al. 1990, 
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Coulson et al. 1994, Gangui & Perivolaropoulos 1995). Moreover, it has been 
shown (Scherrer & Schaefer 1994) that, for a wide class of models leading to 
non-Gaussian density perturbations, the corresponding temperature fluctu- 
ation field, as induced by the Sachs- Wolfe effect, is nearly Gaussian thanks 
to the Central Limit Theorem. Thus, a detailed analysis of the predictions 
for various models, taking into account all the relevant effects, is necessary 
to find out the best tools to distinguish among them. 

The analysis of the three-point correlations in the COBE DMR two-year 
anisotropy maps shows evidence for a non-vanishing signal in the data but 
at a level consistent with a superposition of instrumental noise and Gaussian 
CMB fluctuations (Hinshaw et al. 1995). As the noise level will diminish 
rapidly with additional data, in the four-year map the sensitivity is expected 
to be limited by the cosmic variance. 

Even for the case of primordial Gaussian curvature fluctuations, the 
non-linear gravitational evolution gives rise to a non-vanishing three-point 
correlation function of the CMB. It has been argued by Luo & Schramm 
(1993) that the amplitude of this effect is several orders of magnitude larger 
than that predicted by Falk, Rangarajan & Srednicki (1993) for a cubic 
self-interacting inflaton model, and with a similar angular dependence. A 
more recent estimate of this effect on the skewness of the CMB by Munshi, 
Souradeep & Starobinsky (1995) finds that the amplitude of the gravitational 
and inflationary non-linearity contributions is comparable for typical infla- 
tion models. We provide here a detailed and more general analysis of the 
possible observational consequences of the non-linear gravitational growth 
of initially Gaussian perturbations on the CMB three-point function. 

At large angular scales the anisotropics in the CMB are given by (e.g. 
Martinez-Gonzalez, Sanz & Silk 1990) 



where the first term represents the well-known Sachs-Wolfe effect (Sachs 

6 Wolfe 1967), while the second one corresponds to the Rees-Sciama, or 
integrated Sachs- Wolfe effect (Rees & Sciama 1968). In the previous formula 

7 denotes a direction in the sky and r\ is the conformal time, with rj and 
j] r the present and recombination times, respectively. The contribution of 
the Rees-Sciama effect to the total anisotropy is small in a flat matter- 
dominated universe because the gravitational potential keeps constant in time 
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within the linear regime. It gives however a non-vanishing contribution when 
the non-linear evolution of perturbations is taken into account. This non- 
linear contribution will generate a non-vanishing three-point function, even 
for primordial Gaussian perturbations in the energy density. Second-order 
perturbation theory gives, for the flat matter-dominated case, $(x, r/) = 
$i(x) + $2(x, r/), where the second-order gravitational potential has Fourier 
transform (e.g. Peebles 1980) 

(2) 

and $1 is the linear theory gravitational potential; a(r/) = (f]/r] ) 2 is the scale 
factor. 

The three-point correlation function for points at three arbitrary angu- 
lar separations a, f3 and 7 is given by the average product of temperature 
fluctuations in all possible three directions with those angular separations 
among them. The general expression is given by Gangui et al. (1994). In 
this paper, for simplicity, we will restrict ourselves to the collapsed case, cor- 
responding to the choice a = f3 and 7 = 0, that is one of the cases analysed 
for the COBE DMR data by Hinshaw et al. (1994, 1995) (the other is the 
equilateral one, a = f3 = 7). The collapsed three-point correlation function 
of the CMB is given by 

C M ^J^J ^AT( 7l )AT 2 (7 2 M(7i • 72 - cos a). (3) 

For a = 0, we recover the well-known expression for the skewness, 63(0). By 
expanding the temperature fluctuations in spherical harmonics AT(#, <p) = 
Him O'lmXtmiP i^p)-, we can write the collapsed three-point function as 

C z( a ) = ^ ^ 1 ( COS «)^ 1 m 1 «4m2«4m3 W ^ W ^W4^S 2mi ' 

h /-2/-3 ,m\ ,m 2 ,m 3 

(4) 

where Y\)i represents the window function of the particular experiment and 
we follow the notation in Gangui et al. (1994), defining 

nz7r = J (5) 
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which have a simple expression in terms of Clebsh-Gordan coefficients. The 
mean angular bispectrum predicted by a given model can be obtained from 



(a, iimi a i2m2 a*i 3m3 ) = J d^dtt^dtt^Y^^^Y^ 2 *^^ 3 (73) 

(AT( 7l )AT( 72 )AT( 73 )) (6) 

[for a general discussion of the properties of the angular bispectrum see Luo 
(1994)]. The contribution from the Sachs- Wolfe term has been computed 
by Falk et al. (1993), Gangui et al. (1994) and Gangui (1994), accounting 
for primordial non-linearities due to inflaton self-interactions. The leading 
contribution coming from the Rees-Sciama effect is obtained from 



(AT( 7l )AT( 72 )AT( 73 )) = 

2 f^o Q 

^(^(^o^i^o) ^ dri— $ 2 (x, ?? ; 



) 

X=73(»?0-»?) 



+ (71 <-> 73) + (72 <-> 73), (7) 

with T = 2.726 ± 0.01 K the mean temperature of the CMB (Mather et al. 
1994). At the same order in perturbation theory there are also non- vanishing 
contributions of the type ^(^(^i^o^i^^cO^^s^o, Vr))', compared to the 
non-local Rees-Sciama terms, these terms are however suppressed by about 
two orders of magnitude, because of both the ^-dependence of $ 2 and the 
different numerical factors. Using the analytical expression for $ 2 from Eq. 
(2), a straightforward computation leads to 



(AT( 7l )AT( 72 )AT( 73 )) = 

/"fin 



T 3 r d 3 k t r d 3 k 2 po 
189 7 (2tt) 3 J (2tt) 3 



1]t 
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^2 



\2 



^(ki+e 2 )+7(k 1+ k 2 y - (ki+k2)2 (((k 1 +k 2 )-k 1 ) 2 +((k 1 +k 2 )-k 

P$(&i)P$(&2) + (71 <-> 73) + (72 ^ 73), (8) 

with P$(k) the gravitational potential power spectrum. Using this expression 
and well-known integral relations for spherical harmonics, the three angular 
integrations in Eq. (6) can be performed and we obtain for the bispectrum 

4T 13 r t i"no 

T]r il ,32 ,ni ,n 2 
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\P(h\P(h\ , 20fc 1 2 fc 2 2 + 8(k 1 -k 2 ) 2 + 14(fc 1 2 + fc 2 2 )k 1 -k, 
Y n 2 (k 2 )P$(A; 1 )P$(A; 2 ) ( fc 2 + fc | + 2 k x • k ' ' ^ ' 
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where ji are spherical Bessel functions of order L Now it is possible to make 
the integrations in dfty. and perform the summation over m 8 - in Eq. (4), 
using relations among Clebsh-Gordan coefficients 1 to obtain for the collapsed 
three-point function 



C ^ a ) = - ^Jf 6 E (24 + l)(24 + l)W^W^W 4 (440 0|40) 2 P^(cosa) 
/ dwiw^~ 3 / dw 2 W2~ 3 / ok ^ + i (tfi)^ + i (tfi^) 

JO JO Jrir/vo Z 2 2 



W>?0 

.(10) 



J e2+ i(w 2 )J e+ i(w 2 z) 10K + ^ 2 + 5^ *Mn ^ 



One can readily see from the last expression the advantages of our approach: 
any primordial spectrum may be studied, as we have P$(k) = A(kr] ) n ~ 4 . We 
normalize its amplitude to the Q rm s-PS value determined by COBE through 
A^l = (36/5)2 1 -M 7 rg 2 ms _ P5 r - 2 r 2 (2 - n/2)r(9/2 - n/2)/(r(3 - n)r(3/2 + 
n/2)). Moreover, the expansion in multipoles allows us to specify any chosen 
window function Y\)i as well as to subtract the desired multipole contribu- 
tions from our expression, and thus match the different settings of various 
experiments (provided that we stay far enough from the Doppler peak). The 
integrals in Eq. (10) can be numerically evaluated to obtain the predic- 
tion for the collapsed three-point correlation function produced through the 
Rees-Sciama effect, to be measured by the particular experiment. 

We turn now to the computation of the cosmic variance associated to the 
collapsed three-point correlation function. If the fluctuations in the CMB 
temperature were Gaussian, the mean value of the three-point correlation 
function over the ensemble of observers would be zero. However, as we are 
able to perform observations in just one particular sky, this prediction comes 
with a theoretical error bar, or "cosmic variance", that indicates the typi- 
cal values expected for particular realizations of the Gaussian process. Only 
non-vanishing values larger than this can be interpreted enal of intrin- 

sic non-random phases in the distribution of CMB temperature fluctuations. 



1 In particular, we use the identity X2m=-^(^ m 0\£ m) = (21 + l)(5fco- 
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The expected amplitude of the cosmic variance (Cf (a)) at a particular an- 
gular scale can be computed from Eq. (3) 



(c|(q)) s j^j^jX£j^ l( ^_ ma)i ^. t _ ma) 
(Ar(7,)Ar 2 (7 2 )Ar(^)Ar 2 (-, 4 )). (11) 

We assume that the multipole coefficients are Gaussian distributed random 
variables with angular spectrum (a 7 ^ 1 a^*) = Si 1 i 2 S mim2 (iir /5)Ql ms _p S Ci 1 , 
where 

= r(£ + n/2-l/2)r(9/2-n/2) 
1 r(£ + 5/2-n/2)r(3/2 + n/2)' 1 ' 

as it follows from P$(k) oc k n . Using standard combinatorial properties, we 
can compute {Cl{a))a auss from Eq. (11) 

2(9 6 

{ct(a)) Gauss = 19 7"'-y 2 E c ei c e AsK^iM( 2 ^ + + 1) 

P tl (cosa)(P tl (cosa) + P^(cosa) + P 4 (cos a)) (440 0|40) 2 . (13) 

The last expression takes into account only the effect of the cosmic variance 
on the collapsed three-point function; one can easily modify it to allow for 
the experimental noise, by adding its contribution to the angular spectrum. 

We can now apply our formalism to the COBE DMR measurements. 
We consider a Gaussian window function, Wi ~ exp( — \t{t + l)c 2 ), with 
dispersion of the antenna-beam profile a = 3°. 2 (Wright et al. 1992). Figure 
1 shows the expected collapsed three-point function arising from the Rees- 
Sciama effect, for three different values of the primordial spectral index n = 
0.7, 1, 1.3. Note that all three types of perturbation spectra can be easily 
generated in the frame of inflation models (e.g. Mollerach, Matarrese & 
Lucchin 1994, and references therein). In all cases the perturbation amplitude 
has been normalized to two-year COBE DMR data, using the £ = 9 multipole 
amplitude, ai = 9 //K, according to the procedure proposed by Gorski et al. 
(1994), which leads to Q rms - PS = 24.2, 19.5, 15.9 //K, for n = 0.7, 1, 1.3, 
respectively. The top panel is obtained by subtracting from the map only the 
dipole contribution. In the central panel also the quadrupole and the octopole 
have been removed. In the bottom panel all the multipoles up to £ = 9 have 
been subtracted. This procedure allows easy comparison with the recent 
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analysis of the two-year data by the COBE team (Hinshaw et al. 1995). 
Figure 2 contains the corresponding plots for the rras collapsed three-point 
function obtained from our initially Gaussian fluctuation field, {Cl(a))]i^ uss . 
So far we have assumed full-sky coverage; the effect of partial coverage, 
due to the cut at Galactic latitude \b\ > 20°, increases the cosmic variance 
and can be approximately taken into account by multiplying it by a factor 
1.56 (Hinshaw et al. 1994). An analytic estimate (Srednicki 1993, Scott, 
Srednicki & White 1994) gives a somewhat smaller value, 1.23. Note that 
the subtraction of low order multipoles leads to a strong decrease both in the 
signal and in the cosmic variance. For all considered multipole subtractions 
and angular separations, the expected signal stays typically three orders of 
magnitude below the cosmic variance, which makes the effect undetectable. 
The results are found to be weakly dependent on the considered spectral 
index. For comparison, the amplitude of the collapsed three-point function 
produced by non-linearities in the inflaton dynamics (see Gangui et al. 1994) 
is about a factor of ten smaller than that arising from the Rees-Sciama effect. 
We can also evaluate the skewness parameter S = ((AT /T) 3 ) / ((AT /T) 2 ) 2 , 
which has the advantage of being normalization (i.e. Q rms -ps) independent 
for the considered effect; in the case, e.g., of only dipole subtraction and 
n = 1, S ~ —5. Note that, in the same case, the heuristic estimate by Luo & 
Schramm (1993) leads to S ~ 18, while the method by Munshi et al. (1995) 
yields a slight underestimate, S ~ —1. The corresponding rras skewness 
parameter, arising from the cosmic variance of a scale-invariant Gaussian 
field, is S rms 1.3 x 10 4 {Q rms -p S /19.5 //K) -1 . 

One may try to increase the amplitude of the signal-to-cosmic variance 
ratio, by referring to experimental settings, such as FIRS (e.g. Ganga et al. 
1993) and Tenerife (e.g. Hancock et al. 1994), which are more sensitive to 
intermediate angular scales. However, the little increase of this ratio is not 
enough to make it detectable; furthermore, one has to bear in mind that the 
Tenerife measurement is further affected by a large sample variance. 

Finally, let us stress that our calculations apply to measurements of CMB 
anisotropics on scales large enough to be unaffected by the Doppler peak. The 
possible extension of this effect to smaller scales is discussed by Munshi et 
al. (1995), who however argue that also on those scales the cosmic variance 
largely overcomes the intrinsic non-Gaussian signal. 
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Figure captions 

Figure 1 The collapsed three-point function (73(a), in fiK 3 units, as pre- 
dicted by the Rees-Sciama effect, vs. the angular scale a. The dot-dashed 
line curves refer to the n = 0.7 spectral index, the solid line ones to n = 1 
and the dashed line ones to n = 1.3. The top, middle and bottom panels 
represent the cases l m in = 2, 4, 10 respectively. 

l/2 

Figure 2 The rms collapsed three-point function (C^(a))Q auss7 in 10 4 fiK 3 
units, as arising from the cosmic variance of Gaussian fluctuations. Symbols 
are as in Fig. 1. Note that the top, middle and bottom panels, representing 
the cases l m in = 2, 4, 10 respectively, have different vertical scales. 
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